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We report the first direct observation of the strange b baryon (Sj ). We reconstruct the decay 
E^ — > J/ipE~ , with J/ip — » and E~ — > Att~ — > pir~TT~ in pp collisions at yfs = 1.96 TeV. 

Using 1.3 fb _1 of data collected by the DO detector, we observe 15.2 ± 4.4 (stat.)lj'4 (syst.) E^ 
candidates at a mass of 5.774 ± 0.011 (stat.) ± 0.015 (syst.) GeV. The significance of the observed 
signal is 5.5a, equivalent to a probability of 3.3 x 10 -8 of it arising from a background fluctuation. 
Normalizing to the decay Af, — » J/ip A, we measure the relative rate 

o-(S-) x B(E~ -> J/ipE~) _ n 00 _ Lnnn/ „.„. ^+0.09 



<t(A„) x B(A b J/i>A) 

PACS numbers: 14.20.-c, 14.20.Mr, 14.65.Fy 



= 0.28 ± 0.09 (stat.)Ia.M (syst.). 



The quark model of hadrons [l| predicts the existence 
of a number of baryons containing b quarks, with a hi- 
erarchical structure similar to that of charmed baryons. 
Despite significant progress in studying b hadrons over 
the last decade, only the A& (udb) b baryon has been di- 
rectly observed. The (dsb) (charge conjugate states 
are assumed throughout this Letter) is a strange b baryon 
made of valence quarks from all three known generations 
of fermions and is expected to decay through the weak 
interaction. Theoretical calculations of heavy quark ef- 
fective theory Q and nonrelativistic QCD Q predict the 
mass in the range 5.7 — 5.8 GeV [4j. 

Experiments at the CERN LEP e + e _ collider have re- 
ported indirect evidence of the baryon based on an 
excess of same-sign 5~ ' t~ events in jets Interpreting 
the excess as the semi-inclusive — » Z~£~DjX decay, 



the average lifetime of the S fc is 1.42+q'24 P s I n this 
Letter, we report the first direct observation of the 
baryon, fully reconstructed in an exclusive decay. We 
observe the decay — > J/ipE~, with J/ip — > /i + ^t~, 
5~ — > A7r~, and A — * pir~ . The analysis is based on a 
data sample of 1.3 fb _1 integrated luminosity collected 
in pp collisions at >/s = 1.96 TcV with the DO detector 
at the Fermilab Tevatron collider during 2002 — 2006. 

The DO detector is described in detail elsewhere 
The components most relevant to this analysis are the 
central tracking system and the muon spectrometer. The 
central tracking system consists of a silicon microstrip 
tracker (SMT) and a central fiber tracker (CFT) that 
are surrounded by a 2 T superconducting solenoid. The 
SMT is optimized for tracking and vertexing for the pseu- 
dorapidity region \rj\ < 3 (rj = — ln[tan(6>/2)] and 9 is the 
polar angle) while the CFT has coverage for \rj\ < 2. 
Liquid-argon and uranium calorimeters in a central and 
two end-cap cryostats cover the pseudorapidity region 
\rj\ < 4.2. The muon spectrometer is located outside the 
calorimeter and covers the pseudorapidity region |7y| < 2. 
It comprises a layer of drift tubes and scintillator trigger 
counters in front of 1.8 T iron toroids followed by two 
similar layers behind the toroids. 

The topology of E// — > J/i/jE~ — > J/ipAir~ decay (see 
Fig. [T]) is similar to that of the A;, — » J/ifj A decay; there- 
fore, the reconstruction of the J/ip and A and their se- 
lection discussed below are guided by the strategies ap- 



plied to the A;, lifetime measurement in DO [8|. They 
are then validated with simulated Monte Carlo (MC) 
events. The pythia MC program Q is used to g enerate 
signal events while the EvtGen program [lCJ is used 



to simulate S h decays. 



The Ht mass and lifetime are 



set to be 5.840 GeV and 1.33 ps respectively, their de- 
fault values in these programs. The generated events are 
subjected to the same reconstruction and selection pro- 
grams as the data after passing through the DO detector 
simulation based on the GEANT package ll| . MC events 
are reweighted using the weights determined by matching 
transverse momentum (pt) distributions of J/ip, proton 
and pion from the A& — * J/iph — * J/ippn^ decays in 
MC to those observed in the data. 




FIG. 1: Schematic of the — > J/tpE" — ► J/ip An^ — > 
()i + j-L~) (pn~)iv~ decay topology. The A and E~ baryons 
have decay lengths of the order of cm; the has an ex- 
pected decay length of the order of mm. 



J/ip —>■ ji + decays are reconstructed from two oppo- 
sitely charged muons that have a common vertex. Muons 
are identified by matching tracks reconstructed in the 
central tracking system with either track segments in 
the muon spectrometer or calorimeter energies consis- 
tent with the muon trajectory. They are required to have 
Pt > 1.5 GeV and at least one of them must be recon- 
structed in each of the three muon drift tube layers. The 
dimuon invariant mass M(p + pr) is required to be in the 
range 2.5 — 3.6 GeV. In addition, events must have at 
least one reconstructed primary vertex of the pp inter- 
action. If two or more vertices are reconstructed, the 
one closest to the reconstructed vertex (see below) 
is used. Events containing a J/ip candidate are repro- 
cessed with a version of the track reconstruction algo- 
rithm that improves the efficiency for tracks with low px 
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and high impact parameters. Consequently, the efficien- 
cies for Kg, A, and S~ reconstruction are significantly 
increased. Figure Eta) shows the invariant mass distri- 
butions of the reconstructed E~ candidates (see below) 
before and after the reprocessing. The reprocessing in- 
creases the S~ yield by approximately a factor of 5.5. 
For further analysis, J/ip — ► p + ^ candidates are re- 
quired to have mass 2.80 < M{p, + < 3.35 GeV and 
Pt > 5 GeV. The mass windows here and below are cho- 
sen to be approximately ±5<r and the pr requirement 
ensures that the selected J/ip candidates are above the 
sharp turn-on of the detector and trigger acceptances. 

A — > p7r~ candidates are formed from two oppositely 
charged tracks that originate from a common vertex. The 
track with the higher pt is assumed to be the proton. 
MC studies show that this assignment gives nearly 100% 
correct combination. The invariant mass of the pix^ pair 
must have a mass between 1.105 and 1.125 GeV. The two 
tracks are required to have a total of no more than two 
hits in the tracking detector before the reconstructed pir~ 
vertex. Furthermore, the impact parameter significance 
(the impact parameter with respect to the event vertex 
divided by its uncertainty) must exceed three for both 
tracks and exceed four for at least one of them. These 
selection cuts are the same as those in Ref. [8j|. 
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FIG. 2: Invariant mass distributions of the Air pair before the 
reconstruction for (a) the right-sign Air~ combinations be- 
fore and after reprocessing and (b) the right-sign Air~ and the 
wrong-sign A-7r + combinations after reprocessing. The repro- 
cessing significantly increases the H~ yield. Fits to the post- 
reprocessing distributions of the right-sign combination with 
a Gaussian signal and a first-order polynomial background 
yield 603 ± 34 S~'s and 548 ± 31 H + 's. 

The A candidates are then combined with negatively 
charged tracks (assumed to be pions) to form S~ — > A7r~ 
decay candidates. The pion must have an impact param- 
eter significance greater than three. The A and the pion 
are required to have a common vertex. For both A and 
S~ candidates, the distance between the event vertex 
and its decay vertex is required to exceed four times its 
uncertainty. Moreover, the uncertainty of the distance 
between the production vertex and its decay vertex (de- 



cay length) in the transverse plane (the plane perpendic- 
ular to the beam direction) must be less than 0.5 cm. 
These two requirements reduce combinatoric and track 
mismeasurement backgrounds. 

The two pions from S~ — > A-7r~ — > (p7r~)7r~ decays 
(right-sign) have the same charge. Consequently, the 
combination Att + (wrong-sign) events form an ideal con- 
trol sample for background studies. Figure^b) compares 
mass distributions of the right-sign A-7T - and the wrong- 
sign A-7r + combinations. The S~ mass peak is evident in 
the distribution of the right-sign events. A A7r~ pair is 
considered to be a candidate if its mass is within the 
range 1.305 < M(A7r") < 1.340 GeV. 

— > J/ip S~ decay candidates are formed from J ftp 
and S~ pairs that originate from a common vertex and 
have an opening angle in the transverse plane less than 
7r/2 rad. The uncertainty of the proper decay length 
of the J/i/'S - vertex must be less than 0.05 cm in the 
transverse plane. A total of 2308 events remains after 
this preselection. The wrong-sign events are subjected to 
the same preselection as the right-sign events. A total of 
1124 wrong-sign events is selected as the control sample. 

Several distinctive features of the — > J/tp S~ — > 
J/ipAir~ — > (p + fi~) (p7r~)7r~ decay are utilized to fur- 
ther suppress backgrounds. The wrong-sign background 
events from the data and MC signal events are used 
for studying additional event selection criteria. Protons 
and pions from the S - decays of the events are ex- 
pected to have higher momenta than those from most 
of the background processes. Therefore, protons are re- 
quired to have pt > 0.7 GeV. Similarly, minimum px 
requirements of 0.3 and 0.2 GeV are imposed on pions 
from A and S - decays, respectively. These requirements 
remove 91.6% of the wrong-sign background events while 
keeping 68.7% of the MC signal events. Backgrounds 
from combinatorics and other b hadrons are reduced 
by using topological decay information. Contamination 



from decays such as B 



JfipK*- -> J/ipK%ir- and 



B° — > J/iIjK*~-k + — ► J/ip (K s it~)it + are suppressed by 
requiring the S~ candidates to have decay lengths greater 
than 0.5 cm and cos(# co i) > 0.99, as the S~ baryons 
in MC have an average decay length of 4.8 cm. Here 
8 co \ is the angle between the S - direction and the direc- 
tion from the S~ production vertex to its decay vertex 
in the transverse plane. These two requirements on the 
S~ reduce the background by an additional 56.4%, while 
removing only 1.7% of the MC signal events. The con- 
tribution from the fl^ baryon is estimated to be negligi- 
ble. Finally, baryons are expected to have a sizable 
lifetime. To reduce prompt backgrounds, the transverse 
proper decay length significance of the candidates is 
required to be greater than two. This final criterion re- 
tains 83.1% of the MC signal events but only 43.9% of 
the remaining background events. 

In the data, 51 events with the candidate mass be- 
tween 5.2 and 7.0 GeV pass all selection criteria. The 
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mass range is chosen to be wide enough to encompass 
masses of all known b hadrons as well as the predicted 

+ , 



mass of the 3 b baryon. The candidate mass 
is calculated as M(S~) = M(J/tpE~ 
M(Aw~) + M PDG (J/ip) + M PDG (Z-) to improve the 
resolution. Here M(J/^S - ), M(fi + fj,-), and M(An~) 
arc the reconstructed masses while Mpdg^/VO an d 
Mpdg(E~) are taken from Ref. The distribution of 
M(H^) is shown in Fig. [3]^a). A mass peak near 5.8 GeV 
is apparent. A number of cross checks are performed 
to ensure the observed peak is not due to artifacts of 
the analysis: (1) The J/ifjAir + mass distribution of the 
wrong-sign events, shown in Fig.[3lb), is consistent with 
a flat background. (2) The event selection is applied 
to the sideband events of the S _ mass peak, requiring 
1.28 < M(An~) < 1.36 GeV but excluding the 3" mass 
window. Similarly, the selection is applied to the J /if) 
sideband events with 2.5 < M(fi + fi ) < 2.7 GeV. The 
high-mass sideband is not considered due to potential 
contamination from -0' events. As shown in Fig. [3jc- 
d), no evidence of a mass peak is present for cither 
(p + fi~) (jm~)Tr~ distribution. (3) The possibility of a 
fake signal due to the residual b hadron background is 
investigated by applying the final 3^ selection to high 
statistics MC samples of B~ -> J/ipK*~ -> J/ipK^Tt~, 
B" -> J/tpK%, and A b -> J/i/jA. No indication of a mass 
peak is observed in the reconstructed J/ip 3 _ mass dis- 
tributions. (4) The mass distributions of J/ip, H~, and 
A are investigated by relaxing the mass requirements on 
these particles one at a time for events both in the 3^ 
signal region and the sidebands. The numbers of these 
particles determined by fitting their respective mass dis- 
tribution are fully consistent with the quoted numbers of 
signal events plus background contributions. (5) The ro- 
bustness of the observed mass peak is tested by varying 
selection criteria within reasonable ranges. All studies 
confirm the existence of the peak at the same mass. 

Interpreting the peak as 5^ production, candidate 
masses are fitted with the hypothesis of a signal plus 
background model using an unbinned likelihood method. 
The signal and background shapes are assumed to be 
Gaussian and flat, respectively. The fit results in a 3^ 
mass of 5.774 ± 0.011 GeV with a width of 0.037 ± 
0.008 GeV and a yield of 15.2 ± 4.4 events. Unless 
specified, all uncertainties are statistical. Following the 
same procedure, a fit to the MC events yields a 
mass of 5.839 ± 0.003 GeV, in good agreement with the 
5.840 GeV input mass. The fitted width of the MC mass 
distribution is 0.035 ± 0.002 GeV, consistent with the 
0.037 GeV obtained from the data. Since the intrin- 
sic decay width of the 3^ baryon in the MC is negli- 
gible, the width of the mass distribution is thus domi- 
nated by the detector resolution. To assess the signif- 
icance of the signal, the likelihood, C s +b, of the signal 
plus background fit above is first determined. The fit is 
then repeated using the background-only model, and a 
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FIG. 3: (a) The M(H^) distribution of the HJJ~ candidates 
after all selection criteria. The dotted curve is an un- 
binned likelihood fit to the model of a constant background 
plus a Gaussian signal. The (fi + /j,~)Air mass distributions 
for (b) the wrong-sign background, (c) the H~ sideband, 
and (d) the J/tp sideband events. The mass M(J/ip\Tv) — 
M(fi + n~) + Mp OG {J/i>) is plotted for (b) and (c) while the 
mass M(/i + ^~H~)-A/(A7r _ ) + A/pDG(H _ ) is plotted for (d). 



new likelihood L b is found. The logarithmic likelihood 
ratio -y/2 ln(/3 s+ f,/>C;,) indicates a statistical significance 
of 5.5cr, corresponding to a probability of 3.3 x 10 -8 
from background fluctuation for observing a signal that 
is equal to or more significant than what is seen in the 
data. Including systematic effects from the mass range, 
signal and background models, and the track momen- 
tum scale results in a minimum significance of 5.3a and 
a 5~ yield of 15.2 ± 4.4 (stat.)±^ (syst.). The signifi- 
cance can also be estimated from the numbers of can- 
didate events and estimated background events. In the 
mass region of 2.5 times the fitted width centered on 
the fitted mass, 19 candidate events (8 J/-03 - and 11 
J/tp 3 ) are observed while 14.8 ± 4.3 (stat.)^;^ (syst.) 
signal and 3.6 ± 0.6 (stat.)^!^ (syst.) background events 
are estimated from the fit. The probability of back- 
grounds fluctuating to 19 or more events is 2.2 x 10~ 7 , 
equivalent to a Gaussian significance of 5.2er. 

Figure 3] shows distributions of the proper decay length 
for the 19 candidate events, the 3^ sideband events, 
and the MC signal events plus estimated background 
events. The distribution of the candidate events agrees 
well with that expected from the SjJ" signal while the side- 
band events have a lower mean proper decay length. Due 
to the use of lifetime information in the event selection, 
a lifetime measurement is not made in this Letter. 

Potential systematic biases on the measured 5^ mass 
are studied for the event selection, signal and back- 
ground models, and the track momentum scale. Vary- 
ing cut values and using a multivariate technique of 
different variables for event selection leads to a maxi- 
mum change of 0.020 GeV in the 5^ mass. Subtract- 
ing an estimated statistical contribution to the change, 



7 



? 14 

o 12 
w 

£j 10 
o 



c 

> 

LU 



D0, 1.3 fb" 1 



Data signal 

Data sideband 

MC signal + data bkgd 




0.1 0.2 0.3 0.4 0.5 
Proper decay length (cm) 

FIG. 4: The distribution of the proper decay length in the 
transverse plane of the 19 candidate events in the ±2.5<r signal 
mass window along with that of the events in the sidebands, 
defined to be 5a away from the fitted mass. Also shown is the 
expected distribution from 14.8 MC signal events plus 3.6 
background events. The distribution of the sideband events 
is scaled to the number of events in the signal mass window. 
Kolmogorov-Smirnov tests indicate that the distribution of 
the signal events is favored over that of the sideband events 
with respect to the MC expectation by a ratio of five to one. 



a conservative ±0.015 GeV systematic uncertainty is as- 
signed due to the event selection. Using double Gaus- 
sians for the signal model, a first-order polynomial for 
the background model, or fixing the mass resolution to 
that obtained from the MC E^ events all lead to neg- 
ligible changes in the mass. The mass, calculated us- 
ing the world average values [l[ of intermediate particle 
masses above, is found to have a weak dependence on 
the track momentum scale. This has been verified using 
the Ab — > J/tpA and B° — ► J/tpKg events observed in 
the data. A systematic uncertainty of ±0.002 GeV is as- 
signed, corresponding to the mass difference between our 
measurement and the world average [l[ for the Af, and B° 
hadrons. Adding in quadrature, a total systematic uncer- 
tainty of ±0.015 GeV is obtained to yield the measured 
S^mass: 5.774 ± 0.011 (stat.) ± 0.015 (syst.) GeV. 

The E^ a x B relative to that of the A;, baryon is 
calculated using 



a(E-) x B(E- -> JfrE-) _ e(A b - J/^A) ^ 



ct(A{,) x B(Ab — > J/tpA) e (Sj -» J/^3~) N Ab 

where N s - and N\ b are the numbers of E^ and Ah events 
reconstructed in data. Analyzing the same data and us- 
ing the similar event selection criteria and fitting pro- 
cedure as the E^ analysis, a yield of 240 ± 30 (stat.) ± 
12 (syst.) A(, baryons is determined. The efficiencies to 
reconstruct the decays, e(E^) and e(At>), are determined 
by MC simulation, and the efficiency ratio, e(Af ) )/e(S h "), 
is found to be 4.4 ± 1.3. The uncertainty on e(Ab)/e(E' b ) 
arises from MC modeling (27%), MC statistics (10%), 
the reconstruction of the additional pion in the E^ de- 
cay (7%), and the E^ mass difference between data and 
MC (5%). The largest component, MC modeling un- 
certainty, is due to the difference in the efficiency ratio 
with and without MC reweighting. The efficiency ratio 
is found to be insensitive to changes in A;, and E^ pro- 



duction models. Many other systematic uncertainties on 
the efficiencies themselves tend to cancel in the ratio of 
the efficiencies. We find a relative production ratio of 
0.28 ±0.09 (stat.)iaol (syst.). 

In summary, in 1.3 fb -1 of data collected by the DO 
experiment in pp collisions at i/s — 1.96 TeV at the Fer- 
milab Tevatron collider, we have made the first direct 
observation of the strange b baryon E^ with a statistical 
significance of 5.5a. We observe the decay mode E^ — > 
J/ipE~ with J/ip — > n + fi~ , E~ — > An~ — > pn~TT~. 
We measure the E^ mass to be 5.774 ± 0.011 (stat.) ± 
0.015 (syst.) GeV and determine its a x B relative to that 
of the A h to be 0.28 ± 0.09 (stat.)±o.ol (syst.). 
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